# Data release for reproducing IceCube Upgrade neutrino oscillation sensitivities
This data release contains the chi2 maps used to generate expected sensitivities for the atmospheric neutrino oscillation parameters, corresponding to Fiugres 11 and 15 in [arXiv:2509.13066](https://arxiv.org/abs/2509.13066). Please refer to the publication for a detailed explanation of the sample and analysis. All examples in this data release assume the normal ordering. Please note: this is a simulation-only study demonstrating the expected performance of the IceCube Upgrade and does not represent a measurement to real data!

---

## Files
Along with this readme, you'll fine three chi squared maps and one example python script:
* `modchi2map_nufitwoSK.csv`: chi2 map used to produce Figure 11 (left), in which the injected truth point is the best fit point from [NuFit 5.2 w/o SK (upper panel)](http://www.nu-fit.org/sites/default/files/v52.tbl-parameters.pdf), sin^2(theta\_23)=0.572 and delta m^2\_32=2.43×10^{-3} eV^2. 
* `modchi2map_nufitwSK.csv`: chi2 map used to produce Figure 11 (right), in which the injected truth point is the best fit point from [NuFit 5.2 w/ SK (lower panel)](http://www.nu-fit.org/sites/default/files/v52.tbl-parameters.pdf), sin^2(theta\_23)=0.451 and delta m^2\_32=2.43×10^{-3} eV^2.
* `modchi2map_icecube.csv`: chi2 map used to produce Figure 15, in which the injected truth point is the best fit point from the latest published IceCube DeepCore oscillation result [Phys. Rev. Lett. 134, 091801 (2025)](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.134.091801), sin^2(theta\_23)=0.54 and delta m^2\_32=2.40×10^{-3} eV^2.
* `example.py`: example script demonstrating how to load the csv files, plot the chi squared map, and extract the 90% sensitivity contours shown in Figures 11 and 15.

---

## CSV file contents
Each csv file contains the values for sin^2(theta\_23), delta m^2\_31 in units of 10^3 eV^2, and the metric value at each scan point. The sin^2(theta\_23) values tested are `[0.40, 0.42, 0.44, 0.46, 0.48, 0.50, 0.52, 0.54, 0.56, 0.58, 0.60]`, and the delta m^2\_31 values are `[2.40, 2.42, 2.44, 2.46, 2.48, 2.50, 2.52, 2.54, 2.56, 2.58, 2.60]`. Using all combinations of these two parameters, there are a total of 121 scan points in each file. To convert to delta m^2\_32, we use the value of delta m^2\_21=7.42×10^{-5} eV^2 from NuFit 5.2. The metric value is the modified chi2 value defined in the publication in Equation 1.
